Aging under Shear: Structural Relaxation of a Non-Newtonian Fluid 
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The influence of an applied shear fleld on the dynamics of an aging coUoidal suspension has been 
investigated by the dynamic light scattering determination of the density autocorrelation function. 
Though a stationary state is never observed, the slow dynamics crosses between two different non- 
equilibrium regimes as soon as the structural relaxation time Ts approaches the inverse shear rate 
7~^. In the shear dominated regime (at high 7 values) the structural relaxation time is found 
to be strongly sensitive to shear rate (r^ ~ 7~^) while aging proceeds at a very slow rate. The 
effect of shear on the detailed shape of the density autocorrelation function is quantitatively de- 
scribed assuming that the structural relaxation process arises from the heterogeneous superposition 
of many relaxing units each one independently coupled to shear with a parallel composition rule for 
timescales: 1/r -^ 1/r + A'y. 
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One of the most peculiar and intriguing behaviour of 
soft materials is the strong sensitivity of their flow prop- 
erties to the application of an even slight external de- 
formation 1]. When a liquid flows in a steady shear 
state, the inverse shear rate 7"^ introduces a new rele- 
vant timescale in the dynamics. On a microscopic level a 
non-Newtonian character reflects a competition between 
7^^ and the natural timescale of those particle rearrange- 
ments controlling the macroscopic flow properties. For 
moderate shear rates, i.e. those occurring in the macro- 
scopic world, mainly the slow modes are affected by shear 
and give rise to complex rheological behaviour. Under- 
standing the physical mechanism governing the interac- 
tion between slow dynamics and shear has been the sub- 
ject of great theoretical and numerical efforts in recent 
years. 

On the theoretical side, some of the most powerful 
tools for the investigation of slow dynamics in complex 
systemSjSuch as mode coupling theory |2, Q , mean-field 
models yh trap models |5j and molecular dynamics sim- 
ulations la, S I3 J have been recently extended to account 
for the presence of shear. The general picture emerging 
from these studies indicates that the structural dynam- 
ics is very sensitive to even moderate shear rates when- 
ever 7""'^ becomes of the same order of the characteris- 
tic timescale for structural rearrangements. In partic- 
ular, even starting from non equilibrium, aging states, 
the presence of shear ensures the existence of a station- 
ary state whose correlation functions decay to zero on a 
timescale governed by 7"^. 

Despite the growing amount of numerical and theoret- 
ical work investigating the shear-influenced slow dynam- 
ics, an experimental microscopic counterpart is still rel- 
atively poor. Evidences for a shear dependent structural 
relaxation time have been obtained by diffusing wave 
spectroscopy [3, ll3 and Light Scattering Echo experi- 
ments lll|, though it is not straightforward to deduce the 
microscopic dynamics from the observables probed in the 
multiple scattering regime. No attempt has been made 
up to now to use dynamic light scattering (DLS) in the 



single scattering regime, a technique that directly probes 
the density autocorrelation function, which plays a cen- 
tral role in both theoretical and numerical approaches. 

In this Letter we investigate the evolution of the den- 
sity autocorrelation function of an aging colloidal sus- 
pension subject to a steady shear flow. The sam- 
ple is an aqueous suspension of Laponite, a highly 
thixotropic liquid which undergoes structural arrest on a 
timescale which strongly depends on concentration and 
ionic strength [1^ and that can be as long as few months 
[13 ■ We found that the aging dynamics displays two 
different regimes whose boundary is marked by the con- 
dition Ts7 '^1. As long as the characteristic relaxation 
time Ts is small on the time-scale 7~^, aging is unaf- 
fected by the presence of shear. During aging dynamics 
slows down, and when Ts becomes of the order of 7, the 
system enters a shear dominated regime where aging is 
strongly reduced and the structural relaxation time is 
very sensitive to 7. The intermediate scattering func- 
tions, characterising the slow non-equilibrium dynamics 
of the sheared sample, are well described assuming an 
heterogeneous scenario where the complex dynamics re- 
sults from the superposition of relaxing units each one 
independently coupled to shear rate with a parallel com- 
position rule for timescales: 1/r — > 1/r -I- A'^. 

The system consists of an aqueous suspension of 
Laponite RD, a synthetic layered silicate provided by La- 
porte Ltd. Particles are disk shaped with a diameter of 
25 nm and 1 nm thickness. Laponite powder is dispersed 
in ultrapure water at 3% wt concentration, stirred for 
~ 30 min and then filtered through a Millipore Millex- 
AA 0.45 /im filter. The obtained suspension, which is 
optically transparent and initially liquid, is loaded into 
a home made cone and plate shear cell having a fiat op- 
tical window as the static plate. Incident laser beam 
(He-Ne, 10 mW) and scattered light pass through the 
same optical window. The scattered light is collected by 
a mono-mode optical fiber, detected by a photomultiplier 
and analyzed by a home made software correlator, after 
being optionally mixed with a coherent local oscillator 




FIG. 1: Normalized intensity autocorrelation function for ten 
equally spaced waiting times between 0.1 and 15 hours. Top 
frame refers to aging without shear while bottom frame refers 
to aging under shear with 7 = 223 s~^ 



field. The scattering geometry is fixed (scattering vector 
q = 25 /im~^). The possibility of choosing an heterodyne 
correlation scheme enables direct access to the detailed 
velocity profile in the shear cell. Wall slip and distortion 
from linearity in the velocity profile are found to be neg- 
ligible in the whole shear rate range here investigated. 
The non-equilibrium dynamics of Laponite suspensions 
fiowing in a "steady" shear state is then constantly mon- 
itored through the normalized intensity autocorrelation 
function g(2)(i^^i) = (/(q, i^)/(g, t^ -h i))/(/(Q, i„))2. 

Where (..) indicates temporal average over the acqui- 
sition time T which is always much longer than the 
characteristic slow relaxation time r {T/t > 10^) but 
also much shorter than the time one should wait be- 
fore changes in r are significant (T dlogr/dt^ < 10^^). 
In the single scattering regime and within the Gaussian 
approximation g^'^^t^^t) = 1-1- \Fg{tw,t)\'^ 'l4\ where 

Fg{tw,t) ^ {P-q{tw)Pq{tw + t))/{P-q{tw)Pq{tw)) isthe 

intermediate scattering function of the colloidal parti- 
cles. In order to avoid purely geometric decorrelation of 
scattered light, due to both advection and transit time 
effects |l5j , steady shearing is interrupted every ten min- 
utes and the correlation functions are collected in short 
(few tens of seconds) time intervals during which the ro- 
tor is stopped (7 — 0). The robustness of the results 
with respect to changes in acquisition protocol has been 
checked. The top frame of Fig. ^ shows the evolution of 
g^'^\tw,t), in the absence of shear, for an evenly spaced 
set of waiting times t^ spanning the interval 0.1 — 15 



hours. As already observed in [17|, aging dynamics dis- 
plays a two step relaxation scenario: i) a fast exponential 
relaxation process related to single particle diffusion and 
whose characteristic time Tf remains of the same order 
of magnitude found in very dilute suspensions; ii) a slow 
stretched exponential process related to cooperative rear- 
ranging motions and whose characteristic time Tj grows 
exponentially fast with waiting time tw while becoming 
strongly stretched (stretching parameter (i down to 0.2). 
In the bottom frame of Fig. Q] we reported the evolution 
of g^'^\tw,t), observed on the same set of tw reported in 
top frame, when a shear rate 7 = 223 s~^ is applied. 
Though a stationary state is never reached in the ob- 
servation time window, for long enough waiting times, 
the slow relaxation time grows slower than exponentially 
(the correlation functions come closer), while the shape 
of the relaxation function approaches a constant profile 
(constant stretching exponent). The two step decay for 

^ q y^w 1 ^j- 

Fq{t^,t) = /exp [-{t/T,f] +(!-/) exp [-t/rj] 

where all parameters (/, r^, (3, Tf) depend on t^ and on 
q, provides a very good fit for all the correlators corre- 
sponding to the different 7, t^, values here investigated. 
The presence of a shear induced crossover becomes more 
evident if we plot the fitting parameters (t^) (average 
slow relaxation time Ts/ P r(l//3), where F is the Euler 
gamma function) and (3 (stretching exponent) as a func- 
tion of tyj for different shear rates as in Fig. |21 As long 
as the system's dynamics is fast on the timescale 7"^ 
(indicated by arrows in the plot), non-equilibrium dy- 
namics takes place as if shear were not present: both 
the relaxation time Tg and the corresponding stretch- 
ing exponent j3 evolve with waiting time following the 
7 = curve (full circles). The presence of shear starts 
to affect the dynamics as soon as Ts becomes larger than 
7~^: the growth of t^ is dramatically reduced (even if 
not completely stopped) and the stretching parameter 
j3 behaviour flattens. Though aging is never completely 
absent, the slow relaxation dynamics, for a given wait- 
ing time, appears to be very sensitive to 7 being faster 
and less stretched as 7 increases. If we suddenly increase 
7, rejuvenation is observed leading to faster relaxation 
time and smaller stretching (higher (3). Crosses in Fig. 
121 are the average relaxation times and stretching expo- 
nents of two rejuvenated samples obtained by the two 
subsequent shear rate jumps: 63 s~^ — > 223 s~^ and 
223 s-i ->446s-i. 

It is also instructive to look at the values of (rs(7, tw)) 
for constant i^, as a function of 7, as reported in Fig. 
13 The grey region represents the half-plane (tj) < 7^^ 
where slow dynamics takes place on a timescale shorter 
than 7~^. This region is not affected by the presence of 
shear: the fluid is Newtonian and, similarly to viscosity, 
{ts) is not dependent on 7. On the other hand shear 
plays an important role in the complementary half-plane 
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FIG. 2: Average slow relaxation time (ts) and stretching ex- 
ponent /3 as a function of waiting time tw during aging under 
different shear rates 7: (a)446, (o)223, (o)67, (0)22 s"\ 
Solid symbols (•) refer to aging without shear. Arrows in top 
frame indicate the 7"^ values corresponding to each curve 




FIG. 3: Average slow relaxation time (ts) as a function of 
shear rate 7 for seven evenly spaced waiting times between 
and 30 hours. Grey area represents the half-plane (ts) < 7~^. 
Dashed line is a fit to a 7" power law. 



where (ts) displays a strong sensitivity to 7. The non- 
Newtonian behaviour in the upper half-plane of Fig. |2| 
resembles the same ■j~^ power law (dashed line) observed 
in rheological measurements of Laponite viscosity [Sj . 

The whole scenario depicted above provides a micro- 
scopic counterpart of the strong shear-thinning behaviour 
observed in rheological studies of Laponite [2| and of 
mairy other soft materials. In both real [l| and simulated 
1^,13 hquids, the viscosity crossover from a Newtonian to 
a non- Newtonian regime (power law dependence on 7), 
is usually found to be described by scaling laws such as: 



^7(7) 



l+7Tr, 



(1) 



Since viscosity is related to structural relaxation, one 
could think, as suggested in |g , of a dynamical analogue 
of (^ in the form: 
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r(7) t(0) 



Aj 



(2) 



Which has also the advantage of an easily readable mean- 
ing: shear rate provides a parallel relaxation channel 
to the system which becomes predominant as soon as 
1/t(0,T) ^ A'j. This simple relation is indeed found to 
work remarkably well in computer models for supercooled 
liquids 6]. Using the same kind of reasoning, one would 
expect that even in an aging sample, as soon as the un- 
perturbed relaxation time grows large enough, shear rate 
will fix the relevant timescale and the system will be- 
come stationary even in the non ergodic phases. This is 
actually what has been observed in a number of recent 
numerical and theoretical papers, and also by rheological 
studies. 

We are now in the position of testing |j21l by directly 
investigating the complete time behaviour of density dy- 
namics in the presence of a shear flow. At a first glance 
one would conclude that Fig. |21is in contradiction with 
(O and the expectation that shear stops aging. It is, in- 
fact evident that relaxation time continues to grow, even 
if no longer exponentially fast, at least for more than 
one order of magnitude since it first "feels" the shear 
field. On the other hand, as shown in Fig. 12 when the 
crossover occurs, the value of the stretching exponent is 
about 0.5 or less. This means that the unperturbed dy- 
namics occurs on a broad spectrum of timescales (span- 
ning two decades at least) and that {ts) only represents 
an average relaxation time. Many approaches to slow dy- 
namics in complex systems |l6| suggest that such a broad 
spectrum of time scale arises from the heterogeneous 
character of slow dynamics. Assuming this heterogeneous 
scenario, shear is only effective in interrupting aging for 
those timescales growing longer than a fixed, shear de- 
pendent, cutoff. In other words we expect that (0) holds 
separately for every timescale composing structural re- 
laxation. The above scenario qualitatively accounts for 
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FIG. 4: Ratio F^{tw,t)/Fg{tiu,t) as a function of -yt, where 
Fq{tw, t) is the intermediate scattering function measured af- 
ter a waiting time tw during aging with apphed shear rate 7. 
o: 7=446 s-\ t„=10 h U: ^=AAQ s''^ , C = 14 h A; 7=223 
s~^, tu) = 14 h o: 7=67 s~^, iu,=14 /i. Sohd line is an expo- 
nential fit expl-A-yt] with A = 0.22. 



both the observations of a slowly growing (non saturat- 
ing) structural relaxation time after the crossover and a 
narrowing in the distribution of timescales with respect 
to the unperturbed case. We want to push now the above 
observations on a more quantitative ground. Assuming 
that Q holds separately for each timescale and calling 
G{ty^,T) the unperturbed distribution of timescales 

/■oo 

Fqitw.t)^ G(i^,T)exp[-VT]dT (3) 

we would expect that shear affects (O as follows 

/•oo 

F^{t^,t) = / G(t^,r)expH(l/T + A7)]dT (4) 
Jo 

(5) 



F°{t^,t)exp[-tA^] 



In other words if we divide the correlation function mea- 
sured at a given waiting time t^ and shear rate 7 by the 
corresponding unperturbed (7 = 0) correlation measured 
at the same tw and plot the result as a function of jt, we 
should obtain the master curve exp[— A7t]. 

In Fig. 21 we report the result of such a procedure 
obtained for three different values of 7 and two wait- 
ing times. All curves collapse on the same master curve 
which is well represented by a simple exponential with 
A — 0.22 (solid line). Similarly, we could predict the 
shape of relaxation for given i„ and 7 by simply multiply- 
ing the unperturbed correlation function for the same tw 
by the function exp[— ^7*]. If we do this for the 7 values 
here investigated and fit the result with a stretched expo- 
nential we obtain a prediction for Tsitwjj) and I3{tw,j)- 
The results of such a procedure, for those tw values 
where F^{tw,t) is available, are shown in Fig. Oas solid 
lines. The overall agreement with the directly measured 
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FIG. 5: Average relaxation time (top frame) and stretching 
parameter (bottom frame) as a function of waiting time for 
difi'erent applied shear rates. Symbols are data in Fig. ^ solid 
line is the prediction of the model described in the text. 



data points (open symbols) is very satisfactory and sup- 
ports the picture of the slow non-equilibrium dynamics 
of Laponite as an heterogeneous superposition of relaxing 
units each independently coupled to shear. 

In conclusion, we investigated the effect of shear on the 
non-equilibrium structural dynamics of an aging colloidal 
suspension of Laponite. The presence of a shear flow 
strongly affects dynamics as soon as structural relaxation 
enters the timescale 7"^. In this shear dominated region 
the shear rate dependence of the average slow relaxation 
time (ts) is well approximated by the power law 7"" with 
a ^ 1. The effect of shear on the detailed shape of the 
intermediate scattering function can be well described 
assuming that the slow relaxation process arises from the 
heterogeneous superposition of many relaxing units each 
one independently coupled to shear rate with a parallel 
composition rule for timescales: l/r -^ 1/r + Aj. The 
authors wish to thank F. Zamponi for helpful discussions. 
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